Anderson localization in correlated fermionic mixtures 
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A mixture of two fermionic species with different masses is studied in an optical lattice. The heavy 
fermions are subject only to thermal fluctuations, the light fermions also to quantum fluctuations. 
We derive the Ising-like distribution for the heavy atoms and study the localization properties of 
the light fermions numerically by a transfer-matrix method. In a two-dimensional system one- 
parameter scaling of the localization length is found with a transition from delocalized states at 
low temperatures to localized states at high temperature. The critical exponent of the localization 
length is v » 0.88. 

PACS numbers: 03.75.Ss, 67.85.-d, 71.30.+h 



The question of Anderson localization in an ultra cold 
gas has attracted considerable attention recently by a 
number of experimental groups [H, 0, 0] . Although the 
phenomenon itself has been studied in great detail over 
the last 50 years by many theoretical groups for various 
physical systems [j, its experimental observation 

has been difficult. One of the reasons is that Anderson 
localization is an interference effect of waves due to elas- 
tic scattering in a random environment (disorder). Real 
systems, however, experience also substantial inelastic 
scattering (e.g. absorption of electromagnetic waves by 
the scattering atoms, Coulomb interaction in electronic 
systems etc.) This may hamper the direct observation 
of Anderson localization significantly. Another reason is 
that random scattering is difficult to control in a real sys- 
tem. This is important in order to distinguish Anderson 
localization from simple trapping due to local potentials. 
It requires some kind of averaging over an ensemble of 
randomly distributed scatterers. 

Ultra cold gases offer conditions, where most physical 
parameters are controllable. Since the atoms are neutral, 
there is no Coulomb interaction, and at sufficiently high 
dilution the interatomic collisions are negligible. More- 
over, a periodic potential (optical lattice) can be applied 
by counterpropagating laser fields. This enables us to 
control the kinetic properties of the gas atoms by creat- 
ing a specifically designed dispersion. Disorder could be 
created by disturbing the periodicity of the optical lat- 
tice. In practice, however, this is not easy because real 
disorder would require infinitely many laser frequencies. 
A first attempt is to study the superposition of two laser 
fields with "incommensurate" frequencies (i.e. the ratio 
of the two frequencies is an irrational number) An 
alternative is to randomize the laser field by sending it 
through a diffusing plate 

Recent progress in atomic mixtures pj has offered an- 
other possibility to create disorder in an atomic system. 
Mixing of two different atomic species, where one is heav- 
ier than the other, creates a situation where the light 
atoms are scattered by the randomly distributed heavy 
atoms 

Si, USED]. An 

optical lattice is applied in order 



to keep the heavy atoms in quenched positions. Due to 
their higher mass, the heavy atoms behave classically in 
contrast to the light atoms, which can tunnel in the opti- 
cal lattice. A crucial question is what determines the dis- 
tribution of the heavy atoms. The most direct distribu- 
tion is obtained by putting atoms randomly in the optical 
lattice "by hand" , each of them with independent proba- 
bility Q. This case corresponds to uncorrelated disorder. 
Another possibility is to fill the optical lattice with both 
atomic species and consider a repulsive (local) interaction 
between them. Then the two species have to arrange 
each other such that the total atomic system presents 
a grand-canonical ensemble at a given temperature and 
a given lattice filling. In the presence of interparticle 
interaction within each atomic species there is a com- 
plex interplay of interaction and localization effect. This 
makes it difficult to isolate the effect of Anderson local- 
ization. In order to avoid interaction within each species 
we choose spin-polarized fermions in an optical lattice. 
Then only the Pauli principle controls the short-range 
interaction within each species, such that the remaining 
interaction is between the different fermionic species. It 
has been shown that then the light atoms are subject to a 
quenched average with respect to a thermal distribution 
of the heavy atoms, and that the distribution is related to 
an Ising-like model [9|, llCj, ll2| . The latter implies (strong) 
correlations between the heavy atoms. For systems in 
more than one dimension there is a critical temperature 
T c at which the correlation length diverges. This system 
provides several interesting features for studying Ander- 
son localization. Although it is a many-body system, 
the light atoms behave effectively like independent (spin- 
less fermionic) quantum particles in a random potential. 
The correlation of the randomness can be controlled by 
temperature, where the correlation length decreases with 
increasing temperature for temperatures T > T c , or by 
the strength of the inter-species scattering. 

In the following we shall study diffusion and Anderson 
localization in the grand-canonical ensemble of two spin- 
polarized fermionic species in one and two dimensions. 
Motivated by recent experimental study on a dilute BEC 
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in d = 1 [H, we consider a realistic scenario, in which 
an initial state is prepared at the center with a trapping 
potential and then it is released by opening (i.e. switch- 
ing off) the trap. The expansion of the wave function is 
observed within our numerical procedure. 

model: c' (c) are creation (annihilation) operators of 
the light fermionic atoms, f* (/) are the corresponding 
operators of the heavy fermionic atoms. This gives the 
formal mapping 6 Li -» c\,c r and 40 K ( 23 Na, 87 Rb) -» 
/J, f r . The physics of the mixture of atoms is defined by 
the asymmetric Hubbard Hamiltonian 
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The effective interaction within each species is controlled 
by the (repulsive) Pauli principle, whereas the interaction 
strength of different atoms is U. If the / atoms are heavy, 
the related tunneling rate is very small. The limit if = 
is known as the Falicov-Kimball model 0, El 0, [ll • 

A grand-canonical ensemble of fermions at the inverse 
temperature (3 = 1/ksT is defined by the partition func- 
tion 

Z = Tre-P*. 

In the FK limit if = the Hamiltonian of the light 
atoms H c ({n r }) depends only on the real numbers {n r } 
(n r — 0, 1), and it is a quadratic form with respect to the 
c operators of the light atoms: 
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+ V(C/n r -/i c )c, t c r . (2) 



This means that the density fluctuations n r = f}.f r have 
been replaced by classical variables n r = 0,1. Thus 
H c ({n r }) describes non-interacting fermions which are 
scattered by heavy atoms, represented by n r . The trace 
Tr c in the partition function can be evaluated and gives 
a fermion determinant: 



^=E e 



-Tr, 



{n r } 



-I3h a 



(3) 



The right-hand side is a sum over (non-negative) statis- 
tical weights. After normalization we can define 



P{{n r }) = ie^S."-det[l 



(4) 



which gives 2{n r } P(i n r}) = 1- Thus P({n r }) is a prob- 
ability distribution for correlated disorder and describes 



the distribution of the heavy atoms. In the strong- 
coupling regime tl/2U 3> 1 the distribution becomes 
an Ising model with nearest-neighbor coupling. At half- 
filling (i.e. fj,f = jit c = U/2) this reads 0] 



P({S r }) oc exp -/3(PJ2U) ^ S r S r 



(5) 



<r,r'> 



where S r = 2n r — 1. 

localization: A trapped atomic cloud may be in the 
state \i), concentrated at the center of the optical lattice. 
After switching off the trapping potential the dynamics 
of the atomic cloud is described by the evolution equa- 
tion \ty t ) = e~ lHt \i). The optical lattice remains present 
during the evolution of the cloud. The expansion of the 
cloud of light atoms is studied numerically. Depending on 
the temperature of the grand-canonical system, we find a 
spreading of the wave function at low temperatures but 
a localized behavior at high temperatures (cf. Fig. [TJ. 
Apparently, there is a critical regime with some critical 
temperature T' c , which separates the spreading behavior 
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FIG. 1: Real time evolution of a wave packet. At low tem- 
perature (T = 0: left panel) it propagates due to the checker- 
board configuration of the heavy atoms. At larger tempera- 
ture (T = 0.2: right panel), the wave packet is localized due 
to a disordered configuration of heavy atoms. 
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from the localized behavior. 

The local density of particles at site r of the optical 
lattice with respect to the state \SSf t ) reads 



(*t\4cr\9 t ) = (i\e iHt 4cre- iHt \i) 



with the initial state |z) = Cq|0). The state |0) is an equi- 
librium state of the entire system and can be expanded 
in terms of energy eigenfunctions and Boltzmann weights 
at the inverse temperature (3: 
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For the FK model this expression can also be written as 
a quenched average with respect to the distribution of 
heavy particles [12| : 



{N r ) = (gl(t)g r0 (t)) f 

with the single-particle Green's function 



(7) 



Grr'(t) = 
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(...)/ is the average with respect to the statistical weight 
of Eq. or Eq. ([5]). For a given configuration {n r } of 
heavy atoms the Green's function can also be expressed 
by eigenfunctions of the single-particle Hamiltonian h c 
in Eq. ||2J| (h c <t>k — e-k4>k)- The spatial properties of 
these eigenfunctions determine the spreading of the av- 
erage density particle density (A^ r ) through the Green's 
function: 



Gro(t) = ^2 I 



-ie k t 



1 + e ~0e k 



(8) 



The denominator represents the Fermi function, since our 
atoms are fermions. According to the localization theory, 
it can be assumed that \<f>k,r\ ~ e - ' 1 "'/^, where is the 
localization length. 

After a Fourier transformation of the time-dependent 
density in Eq. ([7]), the u> = Fourier component of (N r ) 
reads 



N, 



,-2|r|/€ 



(l +e -/3e fc )2 ( 1 + e -/3e fc0 )2 



(r 



oo 



where £ is the largest localization length and ek the cor- 
responding energy level. Thus the expansion of the wave 
packet on large scales is controlled by £. 

The localization length can be studied under the 
change of length scales by considering a finite optical 
lattice of length L and width M ||. In particular, we 
analyze the change of the localization length with respect 
to the width M. For this purpose, we define the reduced 



(or normalized) localization length as A m = £/M and 
calculate this quantity by means of a numerical transfer- 
matrix approach [l6j]. Am either increases (delocalized 
states) or decreases (localized states) with the width M, 
depending on the system parameters (e.g. the inverse 
temperature (3). There can also be a marginal behavior 
(e.g. for a special value /3' c ), where Am does not change 
with M. The latter indicates the existence of a phase 
transition from localized to delocalized states. A quanti- 
tative description of the behavior near 8t can be based on 
the one parameter scaling hypothesis [a Gil]. This states 
that In Am can be expanded in a vicinity of the critical 
point Q' c as |17| 



In A M = lnA c ± A\f3 - j3' c \M 



X/u 



(9) 



For A > the positive (negative) sign corresponds to 
delocalized (localized) behavior. Exponentiation of this 
equation and using £ = \(3 — P' c \~" gives 



A 



A c exp 



±A 
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(10) 



where g is the scaling function. Our numerical transfer- 
matrix approach allows us to determine the critical point 
P' c and the exponent v, depending on the interspecies 
coupling parameter U. 

results: First we analyze a one-dimensional system. In 
this case heavy atoms are always disordered due to ther- 
mal fluctuations. The reduced localization length Am 
decreases with increasing length of the system (cf. Fig. 
[2]) at any temperature. This indicates that all states 
are localized. On the other hand, the localization length 
decreases monotoneously with temperature, as a conse- 
quence of the increasing disorder. Therefore, at suffi- 
ciently low temperature the localization length can be 
larger than the size of a finite system. This could be 
relevant in experiments, where we have a finite optical 
lattice. 

In two dimensions the behavior is more complex. First 
of all, the heavy atoms can form an ordered state at 
low temperatures and a disordered state at high tem- 
peratures 0, [l(J. As long as T > 0, thermal excita- 
tions in the ordered state lead to correlated fluctuations 
of heavy atoms. There is a second-order phase (Ising) 
transition with a divergent correlation length at the crit- 
ical temperature T c . The corresponding distribution of 
heavy atoms provides a complex random environment for 
the light atoms. Our numerical transfer-matrix approach 
finds a transition from localized states at high temper- 
atures to delocalized states at low temperatures. There 
is a critical temperature T' c , where this transition takes 
place. For instance, at low temperatures and half filling 
(i.e. for fjtf = ji c = U/2), the heavy atoms are arranged 
in a staggered configuration with weak thermal fluctua- 
tions. Using the approximated distribution of Eq. 
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the effective spin-spin coupling t^./2U leads to the criti- 
cal temperature T c oc t^/2U. The result for the reduced 
localization length at U = 9 (measured in units of i c ) is 
shown in Fig. [3l All curves cross at j3' c ps 16.5, indicat- 
ing a localization transition. With these parameters the 
Ising transition is at (3 C ps 15.9. Therefore, the localiza- 
tion transition occurs in the ordered phase of the heavy 
atoms. The one-parameter scaling function of Eq. (fT0|) 
with 



A r w 10.9, A ps 0.09, 



0. 



(11) 



fits the data of the transfer-matrix calculation (cf. Fig. 

ED. 

In conclusion, we have discussed a mixture of two 
fermionic species with different masses in an optical lat- 
tice, using the Falicov-Kimball model. The heavy atoms 
are represented as Ising spins and the light atoms as 
quantum particles. The latter tunnel in a random en- 
vironment which is provided by a correlated distribution 
of heavy atoms. The distribution of the heavy atoms is 
given by an Ising-type model, which undergoes a second- 
order phase transition in d = 2 from staggered order to 
disorder. Depending on the dimensionality (d — 1,2) 
of the atomic system and the physical parameters (e.g. 
temperature or interaction strength), the quantum states 
of the light atoms are either localized or delocalized. All 
states of light atoms in a one-dimensional fermionic mix- 
ture are localized. In a two-dimensional mixture these 
states are localized at high temperatures and delocal- 
ized at low temperatures. Such a system can be real- 
ized experimentally as a mixture with two spin-polarized 
fermionic species. 
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FIG. 3: Reduced localization length of light atoms for d = 2 
and U = 9. There is a critical inverse temperature (3' c ps 16.5, 
where an Anderson transition occurs. 
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FIG. 2: The reduced localization length Am of light atoms 
in d = 1 for interaction strength (7 = 9. Am decreases for 
increasing system size M, which indicates localized states. 



FIG. 4: One parameter scaling for d — 2 and U — 9. The 
lower (upper) branch represents (de-) localized states. The 
data from the transfer- matrix calculation approach the scaling 
function of Eq. JTDJl with A c ps 10.9, A ps 0.09 and v ps 0.88 
(dashed curves). 
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